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The rare decay 77 is searched for in 104 fb^^ of data, corresponding to 111 x lO'' BB pairs, 

collected with the Belle detector at the KEKB asymmetric-energy e'^e~ collider. No evidence for the 
signal is found, and an upper limit of 6.2 x 10~^ at 90% confidence level is set for the corresponding 
branching fraction. 

PACS numbers: 13.20.He, 14.40.Nd 



The channel — > 77 is a rare decay of the B° meson 
that is interesting both experimentally, for its remarkably 
clean signature, and theoretically, as a tool for constrain- 
ing physics beyond the Standard Model (SM). The SM 
prediction for the B'^ — s- 77 branching fraction (BF) is 
around 3 x 10~® 0] , and a possible Feynman diagram con- 
tributing to this channel is shown in Fig. ^ Sizable en- 




FIG. 1: A possible diagram contributing to 5° — > 77 at the 
lowest order in the SM. The exchange of a charged Higgs 
boson instead of the W boson could contribute to this process 
in some extensions of the SM. 

hancements of the BF are predicted in many new physics 
models Q ; a typical contribution arising from non-SM ef- 
fects would follow from the replacement of the W boson 
in Fig. n with another charged particle such as a charged 
Higgs boson. The —> 77 channel is also interesting 
because it allows the study of non-trivial QCD dynamics 
in B decay, via a pure non-hadronic final state. 

Experimental limits on the BF have been set by L3 Q 
and BaBar 0]. The BaBar upper limit of 1.7 x 10"^ at 
90% confidence level (CL), obtained with 19.4 fb^^ of 
data, is the most restrictive existing experimental con- 
straint on this channel. 

The present search for the 77 decay is based 

on a data sample of 104 fb~^, which contains 111 x 10^ 
BB pairs, collected with the Belle detector at the KEKB 



asymmetric-energy e+e~ (3.5 on 8 GeV) collider |5( op- 
erating at the T(45') resonance. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a 4-layer silicon vertex detector, 
a small-cell inner drift chamber a 47-layer central 
drift chamber, an array of aerogel threshold Cerenkov 
counters, a barrel-like arrangement of time-of-flight scin- 
tillation counters, and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic 
field. An iron flux-return located outside of the coil is in- 
strumented to detect mesons and to identify muons. 
The detector is described in detail elsewhere 0- 

The 77 events are characterized in the center-of- 

mass (CM) frame by two back-to-back highly energetic 
photons. Photons are selected from isolated clusters in 
the calorimeter that are not matched to charged tracks. 
We require a shower shape consistent with that of a pho- 
ton: for each cluster, the ratio of the energy deposited in 
the array of the central 3x3 calorimeter cells to that of 
5x5 cells is computed, and clusters with a ratio smaller 
than 0.95 are rejected. 

In the Belle detector, a large background for this chan- 
nel is due to the overlap of a hadronic event with en- 
ergy deposits left in the calorimeter by previous QED 
interactions (mainly Bhabha scattering). Such compos- 
ite events are completely removed using timing informa- 
tion for calorimeter clusters associated with the candi- 
date photons. Only photons that are in time with the 
rest of the event are retained. The efficiency of this se- 
lection on signal events is larger than 99.5%. The clus- 
ter timing information is stored in the raw data, and is 
available in the reduced format used for analyses only for 
data processed after the summer of 2004, thus limiting 
the dataset available for this analysis to 104 fb~^. 

Rejection of tt*^ and 77 mesons is of primary importance 
in a search for a purely radiative rare decay of the B" 
meson. All pairs of photons with energy larger than 50 
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MeV and a maximum absolute value of the difference 
between their invariant mass and the 7r° mass of 15 
MeV/c^ are identified as neutral pions. For 77 mesons, 
the minimum energy requirement is 100 MeV, and the 
invariant mass of the two photons is required to be within 
60 MeV/c^ of the mass All pairs of photons passing 
either the tt'^ or 77 selection are removed from subsequent 
analysis. 

The two highest-energy photons are selected and their 
momenta are added to reconstruct the momentum. 
Two variables are used to separate signal events from 
background: Mbc = V^bLn/c^ " /c^ and AE = 
E*B - ^^beam' ^hcrc E^^^^ is the beam energy and E^ 
and p*g are the reconstructed energy and momentum, 
all variables being evaluated in the CM frame. The sig- 
nal is concentrated near A_E — and Mbc equal to the 
B'^ mass. The signal window is therefore defined as 

- 0.25 GeV < AE < 0.15 GeV 
5.272 GeV/c^ < A/bc < 5.288 GeV/c^ 

corresponding to about two standard deviation intervals 
above and below the central values just mentioned. 

The main background for the B^ —>■ 77 channel is due 
to continuum events, mostly coming from light quark pair 
production and fragmentation {uu, dd, and ss, uds for 
short). Two variables that display quite powerful sep- 
aration between signal and continuum background are 
a Fisher discriminant based on modified Fox- Wolfram 
moments and the B° production angle with respect 
to the beam in the CM frame, cosO*^. These variables 
are combined in a likelihood ratio, LR; signal and back- 
ground distributions used to construct the LR are ex- 
tracted from Monte Carlo (MC) samples. In the con- 
tinuum background, the two particles that are recon- 
structed as photons are more abundantly produced at 
low polar angle (6**, measured in the CM frame), while 
the signal photons have a flat distribution in cos 9* . Se- 
lection requirements on LR (LR > 0.92) and on the co- 
sine of the polar angle of the most energetic photon in 
the event (|cos0*| < 0.65) are optimized by maximizing 
-/Vsig/^/^sig + A^bck, where iV^ig (A'^bck) is the expected 
number of signal (background) events in the signal win- 
dow. The expected numbers of events are computed for 
an integrated luminosity of 104 fb~^ and assuming for 
the signal the BF predicted by the SM and for the back- 
ground the prediction of the continuum MC. The above 
requirements reduce the continuum background in the 
signal window by a factor of 55, while retaining 31% of 
signal events. 

The total selection efficiency for signal, evaluated using 
MC events, is 11.7%. In data, seven events lie in the 
signal window. They are shown in the AE-M\^c plane 
in Fig. 13 where the signal window is represented as a 
solid-border rectangle. 

Exclusive backgrounds coming from rare B decays 
have been studied by means of large MC samples and 
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FIG. 2: A_E versus Mbc for data events, selected as described 
in the text. Seven events are in the signal window (solid 
rectangle). The fit window is shown as a dashed rectangle. 



only two channels have been found to give non-negligible 
contributions within the signal window: — s- tt^tt^ and 



ing fraction, BF(i3° 



•qiT" . Assuming the measured B''^ 



branch- 



ttOttO) = 1.45 ±0.29 X 10-6 [k 



and the existing limit on the B'^ — > 7777" branching frac- 
tion, BF(S° t^tt") < 2.5 X 10-6 at 90% CL 10], 0.09 
B^ 7r°7r° events and less than 0.06 rjir^ events 

at 90% CL are expected. 

A two-dimensional extended unbinned maximum like- 
lihood fit is performed on AE and Mbc to extract the 
signal yield. The probability density functions (PDFs) 
for the signal are extracted from the MC simulation. 
The photon energy resolution in the simulation is cor- 
rected to match the resolution measured in a photon test 
beam The signal PDFs are parametrized with a 

Crystal Ball lineshape function for AE and a double 
Gaussian for Mbc. 

For the continuum background, a linear shape is as- 
sumed for AE, with the slope free to float in the fit, and 
an ARGUS 13] function for Mbc, with the slope param- 
eter also free in the fit. The exclusive backgrounds enter 
the fit with the normalization described above. They are 
parametrized with a Gaussian PDF for AE and a double 
Gaussian for Mbc. 

The fit has four free parameters: two slopes and the 
numbers of events of the continuum background and of 
the signal. It is performed within the AE range be- 
tween —0.7 GeV and 0.6 GeV and with Mbc greater than 
5.2 GeV/c^ . The fit window is shown in Fig.[21as a dashed 
rectangle. The projections of the fit result on AE (with 
Afbc in its signal window) and on A/bc (with AE in its 
signal window) are shown in Fig. |3| as solid curves; the 
continuum background is shown as dashed curves, the 
signal as the dark shaded regions, and the exclusive back- 
grounds as the light shaded regions. 



The signal yield is measured to be iVgig = 1.8 
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FIG. 3: Projections of the fit result on /S.E (with Mhc in its 
signal window) and on Mbc (with AE in its signal window). 
The fitting curve (solid line) is plotted with data (circles with 
error bars, drawn as asymmetric Poisson confidence intervals) 
and the uds background (dashed line) . The filled regions rep- 
resent the signal (dark shading) and the B and 
B° — > tjtt" backgrounds (light shading). 

responding to a limit on the BF of 6.1 x 10^^ at 90% 
CL, obtained by integration of the likelihood curve up to 
90% of its total area, and including only the statistical 
uncertainty. 

Several possible sources of systematic uncertainty are 
considered. Uncertainties are included in the likelihood 
function as additional parameters and then integrated 
over their respective ranges by assuming Gaussian proba- 
bility distributions. The largest contribution is due to the 
modelling of the signal shape, which depends on angular 
and energy resolutions of the calorimeter. Uncertain- 
ties on these quantities, evaluated by studying samples 
of Bhabha and e+e^— > 77 events, have been propagated 
to the parameters of the signal PDFs and to the fit re- 
sult. Other contributions are the uncertainties on the 
photon reconstruction efficiency, on event selection (LR 
and cos 9* requirements, tt" and 77 mesons rejection), on 
the number of BB events, on background shapes, and 
on the normalization of the exclusive backgrounds. The 
separate contributions are summarized in Tabled as un- 
certainties on the signal yield. 

TABLE I: Summary of the main systematic sources, ex- 
pressed as uncertainties on the fit signal yield. 



Source 


Syst. 


unc. on A'sig 


Signal shape 




0.37 


Photon rec. efficiency 




0.09 


LR and cos 6* req. 




0.06 


tt" and 77 vetoes 




0.05 


Number of BB events 




0.03 


Background shape and norm. 




0.02 



Inclusion of systematic uncertainties results in the fol- 
lowing upper limit on the BF: 

BF(B° -> 77) < 6.2 X 10"^ at 90% CL . 

In conclusion, a search for the decay B^ 77 has been 



performed in 104 fb~^ of data with the Belle detector. No 
evidence of a signal has been observed and a new upper 
limit has been set, corresponding to an improvement of 
the previous limit of about a factor of three. 
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